Effective-range description of a Bose gas under strong one- or two-dimensional confinement 
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We point out that theories describing s-wave collisions of bosonic atoms confined in one- or two-dimensional 
geometries can be extended to much tighter confinements than previously thought. This is achieved by replacing 
the scattering length by an energy-dependent scattering length which was already introduced for the calculation 
of energy levels under 3D confinement. This replacement accurately predicts the position of confinement- 
induced resonances in strongly confined geometries. 



Many experiments investigating the properties of cold 
atomic gases and Bose-Einstein condensates are now per- 
formed in tightly confining traps, such as tight optical lat- 
tices, leading to systems of reduced dimensionality [ fl [2] [3] 
[4). There are many uses for such confinements. In spectro- 
scopic measurements, they eliminate unwanted Doppler and 
recoil effects J3] |6|. They can also be used to create tun- 
able analogs of condensed matter systems, and give the pos- 
sibility to investigate remarkable many-body regimes in low 
dimensions such as the Tonks-Girardeau gas Q [TOl . 
The theory of s-wave atomic collisions in strongly con- 
fined systems has been established in Refs. [7| and ifTTI 
for 2D and ID confinement, respectively. Both predict a 
confinement-induced resonance of the effective ID or 2D in- 
teraction strength. These predictions rely on a description 
of the atomic interaction in terms of the scattering length 
only. However, in 3D confined systems, it was shown that 
a more refined description is needed for very tight confine- 
ment lfT2l [T3l . Similarly, in 2D confined systems, numerical 
calculations in Ref. [14] showed that the scattering length 
description of Ref. |7| may be insufficient. In this paper, we 
present an accurate analytical description for scattering in ID 
and 2D geometries based on the findings of Refs. lfT3l IT31 . 

We consider a gas of bosonic atoms in an optical lattice 
and assume that there is little tunelling between the lattice 
cells, so that each cell is independent. The atoms in a cell 
are confined by a trapping potential which will be assumed 
harmonic (which is true near the centre of the cell). Let us 
consider a pair of atoms in such a cell. For a harmonic po- 
tential, the centre-of-mass motion decouples from the rela- 
tive motion and the stationary Schrodinger equation for the 
relative motion wave function ip(r) reads: 



2p, 



V* + U{r) + V(r) 



tp(r) = Eip(f). 



(1) 



Here, r = (x, y, z) is the relative coordinate with separation 
r, p the reduced mass, U (r) the isotropic atom-atom interac- 
tion potential, V(r) the trapping potential, and E the relative 
energy. 

For 2D confinement (tube or wave guide geometry), the 
atoms are strongly confined in the xy directions and (almost) 
free to move in the z direction, therefore we set: 



V(r) = V 2D (r) 
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where p = \\p\\ and p is the projection of f on the xy 
plane. For ID confinement (pancake geometry), the atoms 



are strongly confined in the z direction and (almost) free to 
move in the xy directions: 

V(r) ee V 1D (r) = X -iiu?z 2 

Here, lo is the trapping frequency at the centre of the cell 
and we define a — \Jhj(2puf) as the typical length scale 
associated to the trap in the confined directions. 

Any scattering solution ip{r) of Eq. Q is composed of an 
incident wave and a scattered wave. A plane wave basis can 
be used for the incident wave, and the scattered wave can be 
expressed with the noninteracting Green's function G(f, r 1 ) 
of the system for U(r) = 0. Namely, for 2D confinement, 
one has: 



(2) 

where (j> nm (p) denotes the unit-normalised 2D isotropic 
harmonic oscillator eigenstate of principal quantum num- 
ber n and angular quantum number m, and q nm — 

\JjjiE — (2n + 1 + \m\)/a 2 is the wave number of the in- 
cident plane wave. The Green's function reads: 



Jq vl i\z~z'\ 



For ID confinement, one has: 

iP(r) = Vn {z)e^P~ ( G 1D {?^)U{r')i>{7)d z ? (3) 



where (p n (z) denotes the unit-normalised ID harmonic os- 
cillator eigenstate of vibrational index n, and q n is the 
wave vector of the incident plane wave with norm q n — 

"j^E — (n+ \)jc> 2 . The Green's function reads: 



G 1D (r,r) 



ip„(z)ip„(z / )^H^ ) (q u \p - p\) 



where Ha^ is the first Hankel function. 

We first assume that the interaction potential U (r) has a 
finite range r^. This means that there is a separation r ^> 
beyond which the wave function is essentially a solution of 
Eq. ([!]) with U(r) = 0, that is to say a solution of the nonin- 
teracting problem. This is indeed the case for typical atomic 
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Interatomic separation z along the weak direction (units of (If,) Interatomic separation p along the weak direction (units of (i f} ) 



Figure 1: Cut through the function rtp(r) in a confined geometry for a van der Waals interaction with scattering length a = 3 fie , where 
fie = ( 2 fc? 6 ) 1,/4 is the van der Waals length. We have chosen a = 3 fie as an illustrative value; for comparison, a ~ 0.6 fie for 87 Rb, 
a ~ —10 fia for 133 Cs, and a > A fie for 86 Sr. Left panel: function zip(p = 0, z) as a function of z for 2D confinement (tube) with 
a — 1.95 fie- Right panel: function p"4>(p, z — 0) as a function of p for ID confinement (pancake) with a = 1.18 fie. In both panels, 
the s-wave component of the solution to the free-space scattering problem at energy E = h 2 k 2 /2p is represented as a thin black line. For 
r > fie it has the asymptotic form corresponding to Eq. The thick black line represents the ID or 2D wave function in the noninteracting 
region obtained from Eq. |4]( or It is determined by a matching procedure with the free-space scattering wave function, as explained 
in the text. The previous theories 171 11 II were based on a matching with the solution to free-space scattering at zero energy (dotted line), 
which has the asymptotic form corresponding to Eq. ^ for r > fie. The resulting noninteracting ID and 2D wave functions are represented 
by thick grey lines. They do not connect to the zero-energy wave function for the considered confinements. The lengths fie, a and 1/k are 
indicated by arrows. 



interactions which drop off as — Ce/r 6 van der Waals poten- 
tials. Because of this fast drop-off, the wave function (at 
sufficiently low energy) reaches its noninteracting form for 



> ^/J£, where /3 6 



(2 J± c s \ 1/ \ 



20A y 6 , w.. C i C H6 - y—f-r- j OZQ. The separation r 
is therefore on the order of /3g, which ranges typically from 
2 nm to 5 nm. 

We further assume that the interaction potential U scatters 
only s waves. The scattering of partial waves of arbitrary 
order under cylindrical confinement was treated in detail in 
Ref. ifm . Retaining only s-wave scattering is valid in the ab- 
sence of shape resonances and if the cold-collision condition 
kro < 1 is satisfied, where hk = y/2/.iE/h is the collisional 
momentum (see Appendix). 

Under these two assumptions, the noninteracting form of 
Eqs (2|3l is obtained by first taking the limit r ^> r' of the 
Green's function (since r > ro ^> r& ~ r') and then approx- 
imating the remaining integral over f" by 4>v,/j,(0) (or ip v (0)) 
times a quantity Air A which does not depend on the indices 
v or /i appearing inside the Green's funtion (taking into ac- 
count the dependence on v, p would introduce higher-order 
partial wave scattering llTTll - see Appendix). These two steps 
are formally equivalent to taking the Green's function out of 
the integral in Eqs. (2|3 1 and evaluate it at r' = 0. Note that 
this can be achieved at any r by replacing the potential U (r) 
by a regularised contact interaction lfT8l . For 2D confine- 
ment, one finds |fl9l 



r>r ~ — - 

and for ID confinement Q, 



i'. it 



(4) 



r>r 



(5) 



wher the factors A mn and A n are to be determined. 

The wave function t/j can also be expanded in spherical 
partial waves: 



(6) 



1=0 m=- 



where Yt, m are tne spherical harmonics and 9 the angle be- 
tween r and the z axis, (p the angle between p and the x axis. 
At short separations, the confining potential V(r) is negligi- 
ble, so one can expect that within a certain range r\ related 
to the confinement length a (see Appendix), the wave func- 
tion is close to a solution of Eq. ([!]) with V(r) = 0, that is to 
say, a solution of the free-space scattering problem. Beyond 
ro, the partial waves of this solution reach their noninteract- 
ing form which is known to be a combination of regular and 
irregular spherical Bessel functions. Consistent with the s- 
wave approximation, there is no irregular Bessel function for 
I 5^ 0, ie no scattered partial wave. For the s wave (£ = 0), 
we have: 



sin kr 
kr 



a{k) 



cos kr 



for ro < r < n , 



(7) 

where 77 is a normalisation factor to be determined and 
a(k) = ~ tan 5k /k is the energy-dependent .s-wave scatter- 
ing length introduced in Refs. ifTSl [151 {8k is the usual s- 
wave phase shift, related to the s-wave component of the re- 
actance matrix K s {k) = — tan<$fc). This energy-dependent 
scattering length contains all the effects of the interaction on 
the wave function in the region ro < r < r%, and for any col- 
lisional energy E. For moderately tight traps a ^> tq leading 
to small collisional energies, there is a range of r for which 
Eq. (j7]i simplifies to: 



v[i- a 

r 



(8) 
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where a — lim^o a (k) is the scattering length of the poten- 
tial. However, for very tight lattices, a may be close to ro 
and only Eq. |7]) holds. 

The essence of the method used in Refs. Q and |[T9l is 
to assume that a ^> ro and match the noninteracting expres- 
sions Q and Q, respectively, with the free-space expression 
([8]l in the region ro < r < r\ where they are both valid. 
(In Ref. [7 1, this is implicitly achieved by use of a 3D reg- 
ularised contact interaction). By performing the matching 
procedure up to zeroth order in the asymptotic expansion in 
r/cr near r = 0, they obtain two relations yielding the two 
unknown factors r\ and A nm (or A n ). From that knowledge, 
they deduce the effective ID and 2D interaction strengths in 
the quasi-lD or 2D regime (q n , q nm <C k): 



9id = 



92D 



h 2 



C(l/2) 
V2 




(9) 



(10) 



with B w 0.915 and C(l/2)/\/2 -1.033, where £ is the 
Riemann zeta function. The singularity in these expressions 
as a function a, a, or c/ corresponds to the confinement- 
induced resonance. Note, however, that these analytical for- 
mulae are valid only when a is large with respect to ro. 

We stress here that the method can be extended by match- 
ing the s-wave component of the noninteracting expressions 
Q or <|3j with the more general free-space expression (|7]). It 
is straightforward to see that a matching to zeroth order in 
r/cr results in the same conditions as those of Refs. |P71 [T9"1 , 
but with the scattering length a being replaced by the energy- 
dependent a(k) in all the formulae. In particular, a is replaced 
by a(k) in the expressions for rj, A nm , A n , and Eqs. |9p0|. 
Surprisingly, we find that once the two unknown factors r\ 
and A nm (or A n ) are determined this way, the expression Q 
or ^ match the expression ^ up to second order in r /a, ie 
the next two orders automatically match without the need for 
extra parameters. This was checked numerically both in 2D 
and ID confinement, and we give in the Appendix an ana- 
lytical derivation of this result for 2D confinement, based on 
mathematical assumptions which we checked numerically. 
As a result, the expressions Q or Q and dTl) can be matched 
within less than 1% for a > ro. In Fig.fijwe illustrate the 
matching of the wave functions for a confinement a close to 
ro, in the case of a van der Waals interaction U(r) — — %, 
for which ro « /3e, as stated earlier. The figure shows both 
the s-wave component of the solution to the 3D free-space 
problem (which is highly oscillatory for r < ro, and has 
the asymptotic form <J7J for r > r ) and the noninteract- 
ing wave function along either the z or p directions. For 
better comparison, we should plot only the s-wave compo- 
nent of the noninteracting wave function, but this happens 
to be very slowly converging numerically. The inclusion of 
higher-order partial waves thus induces a second-order dif- 
ference (because f/'2,m(r) oc r 2 for small r), but one can still 
appreciate the very good matching of the two functions even 
at these extreme confinements. In contrast, the two functions 
calculated from the original theory do not match. 

The fact that we can gain two orders in the matching pro- 
cedure by using the free-space expression |7} implies that the 
replacement a — > a(k) greatly improves the accuracy and the 



range of validity of the original theory. Although we do not 
have an estimation of the accuracy on g\ and 52, we expect 
the extended theory to give reasonnably good predictions for 
a > r . This improvement brought by the simple replace- 
ment a — > a(k) not so surprising in view of the results re- 
ported in Refs. H3J [B] ED [221 El, which have shown the 
relevance of energy-dependent scattering lengths for the ac- 
curate calculation of energy levels in 3D confined geome- 
tries. It has been long known in other contexts that con- 
tact interaction models are significantly improved when the 
coupling constants are proportional to the energy-dependent 
scattering length (or reactance matrix) [24| . Similar exten- 
sions of |9]) and ( fT0) i were considered in Refs. E51 l26l [271 
in order to take into account the energy dependence due 
to a scattering resonance at low energy. In Ref. ETIl . a 
renormalised contact interaction was used, leading to the re- 
placement of a by the quantity j-^T(k), where T(k) = 

^0-a(k)(l + ika(k))- 1 is the T-matrix. This complex 
quantity is equivalent to the real a(k) at low energy. Here 
we focus on the energy dependence for strong confinement 
even in the absence of any resonance. 

In the quasi- ID (q nm -C k) or quasi-2D (q n -C k) 
regimes, the collisional energy is set by the confinement 
length a and can be relatively high, although the effective 
collisions in the weak directions are cold. An interesting 
consequence of the extended theory is that if the confine- 
ment is strong enough, the collisional energy can probe the 
energy-dependence of a(k). For a standard contact potential 
lfl8l . a(k) is constant and equal to a. However, a more real- 
istic a(k) has some energy dependence. For instance, in the 
effective-range approximation, a(k) has the resonant form: 



a(k) 



1 



T fc 2 c 



(11) 



where r e is the effective range of the potential |28|. This 
approximation works well for short-range interactions with a 
large scattering length a ^> r e . In the case of van der Waals 
interactions, the effective range r e is a simple function of a 

and /3 6 



3 a 



- 1 



(12) 



where a — 2tt(3q/T(1/4) 2 . More elaborate analytical ex- 
pressions of a(k) valid for any a have been derived for van 
der Waals interactions ll30ll . The interest of Eqs. (Hi and 
( 12 1 is that they give a simple two-parameter description of 
the collisions for a wide range of energies. 

To illustrate these ideas, we calculated the ID interac- 
tion strength gin for a van der Waals interaction consistent 
with the Lennard-Jones parameters of the numerical calcula- 
tion reported in Ref. fl4l . The authors observed a differ- 
ence between their numerical calculation and the analytic 
formula (|9]l where a is taken as the zero-energy scattering 
length. They suggested that this difference comes the fact 
that the confinement-induced resonance in gio results from 
a Feshbach resonance with a trap bound state, whose bind- 
ing energy is not predicted accurately by a pseudopotential 
based only on the scattering length. As a result, the for- 
mula Q does not predict the resonance at the right location. 
However, we show in Fig. [2] that the same formula used in 
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Figure 2: The ID interaction strength (left panel) and the 2D interaction strength (right panel) in dimensionless units as a function of the 
ratio of the 3D scattering length a over the confinement dimension a (a is varied and a is fixed). As in Fig. [I] a — 1.95 (5% in left panel 
and a = 1.18 (3q in the right panel. The dots are obtained by numerically solving the problem of two atoms interacting through a Lennard- 
Jones potential in a 3D cylindrically-symmetric trap; dots in the left panel (tube geometry) were taken from Ref. 1141 and confirmed by us, 
while we calculated the dots in right panel (pancake geometry) using the adaptive grid refinement method of Ref. |20|. The dashed curves 
correspond to the formulae ^ and \W) with the zero-energy scattering length a, and the solid curves corresponds to the same formula? using 
the energy-dependent scattering length a(k) instead of a. Here, a(k) has been calculated in the effective-range approximation, Eq. 
which is valid almost everywhere except for small positive scattering lengths where the approximation causes a spurious resonance. 



conjunction with the replacement a — ► a(k) in the effective 
range approximation reproduces the numerical calculations 
very well. This is because the effective range approxima- 
tion is able to reproduce the binding energy of the last bound 
state accurately. The only region where the effective range 
approximation fails is for small scattering lengths a <C (3&, 
where it predicts a spurious resonance, as visible in Fig. [2] 

We also calculated the 2D interaction strength and 
checked that a similar situation occurs in the pancake con- 
figuration. Using the adaptive grid refinement method of 
Ref. Il20ll . we solved the Schrodinger equation ([TJ for a 
Lennard-Jones interaction U(r) — — §f and a cylin- 
drical harmonic trap. The tight pancake limit is obtained by 
setting the ratio of axial and radial frequencies to 400 (thus 
leading to a spatial aspect ratio of 1/20), and the tight con- 
finement scale is set to a = 1.18/36- The parameter C12 
is adjusted to set the number of bound states supported by 
the interaction and the scattering length. From this calcula- 
tion, we obtained the eigenenergies and then used Eq. (21) of 
Ref. BTI to extract the 2D scattering length. We found that it 
shows very little dependence on the number of bound states, 
which can be as low as 2, saving computational efforts. Us- 
ing Eq. (7) of Ref. ED (or Eq. (15) of Ref. l32l ). we could 
then relate the 2D scattering length to the interaction strength 
<?2D for any qo - we chose a qo given by the zero-point mo- 
mentum in the weak direction. Figure plicompares this nu- 
merical g2D with the analytical formula |9| for the same qo- 
Again, the position of the confinement-induced resonance 
for negative scattering lengths [11] is correctly predicted by 
(|9]l provided the energy-dependent scattering length is used. 
As previously, the effective range approximation works well, 
except for small scattering lengths. These results also sug- 
gest that the observation of the resonance may provide useful 
information about the effective range of the interaction. 

In summary, we have shown that the effective ID or 2D 
interactions of ultracold bosons in strongly confined systems 
are governed by 3D collisions at a relatively high energy 
determined by the confinement. The effect of these high- 
energy collisions can be well described by a single quantity, 



the energy-dependent scattering length, up to extremely tight 
confinements. For van der Waals interactions, this quantity 
itself can be expressed in the effective range approximation 
in terms of the zero-energy scattering length and the van 
der Waals length. This parametrized energy-dependent scat- 
tering length leads to an accurate analytic prediction of the 
confinement-induced resonance both in ID and 2D confine- 
ments. 



Appendix 

In this appendix, we show that for 2D confinement the 
matching procedure is effective up to second order in the ex- 
pansion in r/a near r — 0. Without loss of generality, we 
take the wave function to be symmetric around the z axis (ie 
m = in Eq. Q), so that Eq. ([6]) can be written as 

00 

^(r) = J2M r )i\2t+ l)Pt(cos9) (13) 

where Pe(x) is the Legendre polynomial. The basic approxi- 
mation is to assume that the components tpt are proportional 
to those of the free-space solution. With the assumption that 
only the s-wave component is scattered by the potential, we 
set: 

Mr) = Vo<f>no{0)(jo{kr) + ka{k)y {kr)) (14) 

r>r a 



i>&to(r) = rn4>no{0) h(kr) 

r>r 



(15) 



where ji and yi are the spherical Bessel functions, and rjt 
are factors to be determined. (Note that with this definition, 
X] in Eq. ([7]) is equal to 77o</>no(0)). Since only components of 
even I are coupled to the s-wave component by the trapping 
potential V(r), we can discard odd-£ components (they do 
not play a role in the scattering process). An expansion near 



of Eq. ( 13 1 therefore reads 
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ka ( k ) \ ru ^3 
m -24T ) {kr) 



(16) 



This is to be matched with Eq. (j4j). Using the explicit form 
of the harmonic oscillator solutions cf) n o, it can be written as 



Matching the expressions ( 16 1 and ( 17 1 along the z di 



rection leads to the following relations for each order of the 
expansion: 



^no(0)T„ 



y/2a y/2a 2 



2 2 



(17) 



ri a(k) = A n0 

no 



with 



A n 



T„ [u, v, e] = e 2 L n {u ) cosh(2iV n + e) -^-A[u,v,i 

\j2a 



no 



k 2 a(k) A n Q , 2 2\ 
(k a ) 



A [ u > w ' e l = X! 



1^— 



where the L u is the Laguerre polynomial. For Eq. ( 17 1 a 



partial wave expansion appears not feasible, so we resort to 
evaluating the expressions along the p and z axes. Along the 
z axis (p — 0), we find the following expansion 



A[0,u,c] 



1 



+ C(|,e) + (2e-l> + 2C(-i e)i 



(2e — l) 2 — 



1 



(18) 



Along the p axis (z = 0), we find the following expansion 

A[u,0,e] - 1 +C(|,e) + (2e-l)u 



,(2e-l)((l E )-((-i e) 



(2e - I) 2 



(19) 



In these expansions, £ refers to the Hurwitz zeta function 
which we define as follows 



lim 

JV^o, 



N 



1 



(AT + e+ |)- s+1 



(n + e) 5 



for —1 < s < 1, generalizing the definition given in 
Refs. Ifl4ll33l . The first expansion ( 18 1 was found using the 
counter-term method explained in these references with the 

refined counter term f Q + 2 +e ^^^^-ds. The second ex- 
pansion ( [T9] > was guessed from the expected result ( [To} , and 
checked numerically. Because of the very slow convergence 
of the sum in A[u, 0, e] , we could not check the terms di- 
rectly. Instead, we noted that the Laplace transform of A 
with respect to the argument u is related to the Lerch tran- 
scendent $ 1 34 1 : 



J A[w, 0,e}du ■ 



8 +h U+i'2' £ 



and assumed there is a direct correspondence between the 
Laplace transform of the terms in ( |T9] > and the terms in the 
asymptotic expansion of the Laplace transform near s — > 00. 
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while, matching along the p direction leads to the following 
relations 



r) a(k) 



A, 1 
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(27) 
(28) 



IT 



C(I e) + C(-i e) 



24cr 4 



(29) 



Relations |20|2T 



where we make use of 2e — 1 = — k 

and ( |25p6l > are the same and determine A nQ and the s-wave 
factor: 

1 
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"(fc) ft 1 



(30) 



Relations ( |22] > and (|27ji are the same and are both satisfied 
with the previous determination of A n o and rjo. Relations 
p3] > and ( |28| are consistent and give the same determination 
of the d-wave factor: 



»72 = Vo 



( 3a(k) 
\V2k 2 o 



3 Sl 2 ' — Y 
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^ g 2 Q 
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(31) 



However, neither relation ( |24) > or ( |29| i are satisfied with 
the previous determination of A n0 and rj (because of the 
and I), which means that the free-space ap- 



terms 



proximation (14 15 1 is valid up to order r . The error is 



<Pno(0)\Vo a (k) 1 2^pt r3 - The range of r for which this error is 
negligible determines the range [ro, r{\ where a matching is 
possible. Since the error increases with r, we can define r\ 
as the r for which the ratio between the error and the wave 



function ( 16 1 is equal to a certain tolerance a. In the limit of 
large scattering lengths, one finds 

n = (24a) 
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The range [ro, ri] where the two functions can be matched 
to within the tolerance a exists as long as ri > ro, ie 



a > (24a)" 1 / 4 



'■()• 



(32) 



For instance, for a = 1%, one gets a > 1.43ro. This indi- 
cates that the method works even for a confinement length a 
on the order of the range ro. 

We can also check that this is consistent with our 
assumption that only s waves are scattered. Higher- 
order partial wave scattering arises if we take into ac- 
count the dependence of A n o on v. The first correction 
is A n Q — > A„o + (y + i)B n Q, and leads to the term 
Bnov 2 1 'z 3 along the z axis. This term is to be matched 
with the leading-order term of the scattered d-wave term 
(-X2)7?2</>no(0)y2(fcz)(-5)P 2 (cosO), where K 2 is the d- 
wave reactance matrix element. In the absence of shape res- 



onance, Ki is purely determined by the long-range van der 
Waals behaviour of the interaction, K 2 ~ — ^(fc/^) 4 - see 
Eq. (8) of Ref. ED. This leads to: 



B n0 = -r) 2 



v (fc/? 6 ) 4 
20 (kaf 



e?-wave scattering is negligible if \B n o\ <C |^4 n o|- In the 
quasi- ID regime (q„o "C k and k ~ 1/cr), and using the 
value of 772 pi) , one finds: 



at 



which for a wide range of scattering lengths is consistent 
with the condition ( |32| i and the cold collision condition 
kr < 1 given in the text. 
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